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Abstract 

We calculate a correlation function of the Jordan- Wigner operator in a class of 
free-fermion models formulated on an infinite one-dimensional lattice. We repre- 
sent this function in terms of the determinant of an integrable Fredholm operator, 
convenient for analytic and numerical investigations. By using Wick's theorem, we 
avoid the form-factor summation customarily used in literature for treating similar 
problems. 

1 Introduction and main results 



In this paper we investigate the correlation function 

it fnW ^ — A _ a' 



[i.D 



D{X; n, t) = {TO,{t)O],{0)), n = j- f 
of the Jordan- Wigner operator 

Oj = e'^^, (1.2) 

where 

j 

Afj = ^ ^ Qmi Qm — (^m^m- (l-^) 

m=—oo 

The operator Cm (cm) is the fermion creation (annihilation) operator at the m- 
th lattice cite of an infinite one-dimensional (ID) lattice. These operators obey 
canonical anticommutation relation, {cm,cj^/} = 5mm'- The symbol T in Eq. (jl.ip 
stands for the time ordering. The function (jl.ip is a 27r-periodic function of a real 
parameter A. The time evolution of the operators in Eq. (II. ip is generated by the 
free-fermion Hamiltonian 
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Figure 1: Shown is an example of the single-particle energy as the function of quasi- 
momentum k in the first Brillouin zone, — tt < k < n (solid line). For the ground state 
average standing in Eq. (11. ip the /c-space states between —kp and kp are occupied with 
fermions. This is illustrated with the bold line. 

where is the single-particle energy as a function of the Bloch quasimomentum in 
the first Brillouin zone, — vr < k < tt. For simplicity we assume that the equation 
^fc = has only two solutions k = zizkp. In this case, the ground state average (• • • ) 
in Eq. (If .1 p is taken over the Fermi sea of fermions filling the single particle states 
between —kp and kp : 



An example of satisfying Eq. (jf .5p is shown in Fig. [TJ 

The goal of the paper is to derive a Fredholm determinant representation for the 
function (II. ip . This representation was used in our recent paper [l] to investigate 
correlation functions in the ID Bose-Hubbard model with spin. However, in Ref. [1] 
we neither give the explicit expressions (II.6p - (II.I2p nor explain a way they can be 
obtained. It is important to stress that the function (jl.ip appears not only in the 
Bose-Hubbard model with spin. In the regime of strong interparticle repulsion it 
is possible for a broad class of models carrying spin and charge degrees of freedom 
to factorize spin and charge dynamics in such a way that the charge part is fully 
described by the function (jl.ip . This was demonstrated in Refs. [21 [3l [U EJ El [71 [I] 
by representing the original correlation function as a convolution of the charge 
correlation function (jl.ip with some spin correlation function of the Heisenberg 
model (see, for example, Eq. (3.5) of Ref. 0, Eq. (8) of Ref. 0, Eq. (II) of 
Ref. [l], etc.). The authors O [3l [H (H El [7] do not give Fredholm determinant 
representation for the function (jl.ip . analyzing it by other methods. An approach 
complementary to the one of Refs. [HElElllElElElis used in Refs. There ID 

gas of bosons (and fermions) with spin is considered in a limit of infinitely strong 
short-range interparticle repulsion. The explicit expressions for the correlation 
functions are obtained by summing up matrix elements of the transitions to the 
intermediate states (form- factors summation technique) . The results have the form 
of a Fredholm determinant convolved with some weight function. When compared 
with the predictions of Refs. [H El O [H El El E] these results can be unfolded to the 
form of the convolution of some spin correlation function with the function (jl.ip . 
This way one can indeed get the Fredholm determinant representation for the 
function (jl.ip . However, in our opinion, it is easier to calculate the determinant 




k G [—kp, kp] 
otherwise 



(1.5) 
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directly, as explained in our paper. We note that correlation functions similar 
to the one defined by Eq. (II. ip appear in the analysis of many other important 
ID systems. The number of works on the subject is large and we, unfortunately, 
cannot give the exhaustive list of references here. They include the studies of the 
Tonks-Girardeau gas [IHl [HI [111 US] , impenetrable anyons OUS], XY Heisenberg 
chain [TBI [13 [13 [iHl [20l [21] , non-linear quantum hydrodynamics [22] , and the full 
counting statistics of quantum transport [23] . 

The Fredholm determinant representation for the correlation function (II. ip 
reads 

D{X;n,t) =det{i + k), (1.6) 
where K is a linear integral operator 

{kf){k) = r ^K{k,p)f{p) (1.7) 

J-kp 

with the integral kernel 

^ ^ U{k)lAp) -Uk)U{p) _ 2^ _ £_(A:)£^(p)] (1.8) 

tan 

defined on [—kp, kp] x [—kp, kp]- The functions entering this kernel are defined as 
follows: 



and 



sin A 1 — cosAi?(A;) 

z(k) H T-f 

2 ^ ^ 2 z{k) 



(1.10) 



where 
and 

^(A;) =e^("^-«**). (1.12) 



The linear integral operator with the kernel (II. Sp belongs to the class of integrable 
operators, introduced in Refs. [Ml[2S]) see also Ref. [13] for further discussion. For 
the particular case of A = vr the kernel (jl.Sp appears in Refs. [191 [SQl [H] for the XY 
Heisenberg chain and in Ref. [H] for the impenetrable Bose gas. The integrable 
operator obtained in Refs. [3 [9] coincide with ours up to a finite rank operator. The 
continuum limit of the kernel (11. Sp is encountered in the theory of impenetrable 
anyons [15j. 

In some of the above mentioned works (for example, in Ref. [16], see also Ref. |18] 
for the list of references and detailed discussion) Wick's theorem is used to get the 
determinant representation of correlation functions, in others form-factor summa- 
tion technique is employed. Within the form-factor approach the calculations have 
to be carried out in a finite system, and the thermodynamic limit is applied to the 
end result only. In contrast. Wick's theorem approach makes possible to work in 
the thermodynamic limit from the very beginning. Though less rigorous, the latter 
approach enjoys the advantage of being minimalist. In particular, it avoids con- 
structing the determinant representation for the function (jl.ip in the system with 
a finite number of particles and, therefore, significantly simplifies the presentation 
as compared with the form-factor approach. This motivates the use of the Wick's 
theorem in the present work. 
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Figure 2: Shown are two sets of points: Si, Eq. (12.91) . and ^2, Eq. (12.101) . and the 
enumeration of their elements used in Eq. fl2.13p . 

2 Derivation of the results 

Consider a countable set M of points in space-time. For brevity we shall use the 
notation q{x) = Qj{t), c{x) = Cj{t) and c^{x) = cj(t) for the operators at the point 
X = {j,t) from M. For an arbitrary function 7 : — > C consider the average 

CO 

{Tl[[l + j{x)e{x)]) = l + Y, E n 7(x)ct(x)c(x)), (2.1) 

x€M Q=l MqCM xeMq 

where the inner sum is taken over all distinct subsets Mq of ^A containing Q ele- 
ments. Applying Wick's theorem to the average on the right hand side of Eq. ()2.ip 
we get 

(T JJ ^{x)c\x)c{x)) =deiM{MQ). (2.2) 

x&Mq 

In this expression detM(A^g) denotes the determinant of a Q x Q matrix, which 
for a given enumeration xi, . . . ,xq of the elements of Mq is defined by 

[M{MQ)]a,b = l{Xa)G{Xa-Xk), a,b=l,...,Q. (2.3) 

Here 

G{x) = G{j,t) = -(Tc,(t)4), G{j,0) ^ limG(j, -6) (2.4) 



is the T-ordered single-particle Green's function. Combining Eqs. (|2.4p and ()2.2p 
we find 

(T n [1 + lix)Q{x)]) = det(/ + K). (2.5) 

where if is a linear endomorphism of the linear space of complex- valued 

functions on Ai : 

iKf)ix)= ^ jix)Gix-x')fix'), feriM). (2.6) 

x'eM 

Let us now apply the formula (j2.5p to the correlation function (jl.ip . Using the 
identity = qj following from the anti-commutativity of fermionic fields we get 

e*M=i + ^^., ^ = e^^-l. (2.7) 
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The equation (jl.ip is then written as 

D{X; n, t) = (T J] [1 + ^g{x)] J] [1 + 7*^>(x)]) , (2.8) 

where 

Si = {x: x={n+l-a,t), aen*}, (2.9) 

^2 = {x : X = (1 - a, 0), a G N*} (2.10) 

and 

N* = 1,2,3,... . (2.11) 

The right hand side of Eq. ()2.8p can be expressed as a determinant, Eq. ()2.5p . 
where -ftT is a hnear endomorphism of "^(S*! U ^2) = y{Si) © "^(^2). We represent 
as a 2 X 2 matrix 

^ = ( |J1 I22 ) ' ^""^ : ■^('^a) - -^(5/3), «, /3 = 1, 2. (2.12) 

Using the enumeration of elements of Si and ^2 defined in Eqs. (12. 9p and ()2.10p 
and shown in Fig. [2] we write a matrix element of the representation (j2.12p : 

It is convenient to extend this definition to a, 6 € Z by introducing 

rfn _ / ^a+i,fe+i, (a, 6) G Z+ X Z+ 

^^'^'^ ~ \ 0, (a, 6) G (Z X Z)\(Z+ X Z+) ' ^ ^ ^ 

where 

Z+ = 0,1,2,3, ... Z = 0,±l,±2,±3, .... (2.15) 
Since det(/ + U) = det(/ + K) we get for Eq. 1^ 



D{X;n,t) = det{i + U). (2.16) 
Defining the Fourier transform by 

^e^'^^m, f{k) = e-^'^fij) (2.17) 

and applying it to Eq. (|2.14p we express U as a linear integral operator 

dk 

^-U{k,p)f{p) (2.18) 

with the integral kernel 

U{k,p) = r '^y- e-'i^''-''^e''"'e-'P^u{q), (2.19) 



a,b=0 



where 
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The function G entering Eq. ()2.20p is the Fourier transform of Green's func- 
tion (ITU) : 

G{q,t) = -e{t > 0)[1 - ni.(g)]e-*«'* + 9(1 < 0)n^(g)e-^««*, (2.21) 

where .^g is the single-particle energy, Eq. ()1.4p . and npiq) is the Fermi-Dirac 
function (jl.Sp . Performing the summation in Eq. ()2.19p we obtain 



U{k,p) = j ^vik-qHQ-pHq), (2.22) 

where 

^(^) = 13^1^ (2.23) 

and e is an infinitesimal positive constant. 
Next, we note that 

u{q)=v{q)+w{q), (2.24) 

where 

= ( -7° (,) ) 

and 

wiq)=nHq)[j^^^ (2.26) 

The function z{q) is defined by Eq. (jl.l2p . Using the decomposition (j2.24p we write 
U = V -\-W^ where 

V{k,p) = ^rj{k - q)rj{q - p)v{q) (2.27) 

and 

r da 

Wik,p) = J -^r]{k-q)7]{q-p)w{q). (2.28) 

It follows from Eq. ([2:25]) that = Q. This imphes 

(/ + y)-i = /-t> and det(/ + y) = l. (2.29) 
Using the identities (|2.29p we rewrite Eq. ()2.16p as 

D{X-n,t)=det[i + {i-V)W]. (2.30) 
The integral operator W can be represented as a product 

W = AB, (2.31) 
where the operators A and B possess the following integral kernels 

A{k,q)=r^{k-q)a{q) B{q,p) = r^{q - p)b{q) (2.32) 

and 
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Since det(i + AB) = det(j + BA) we rewrite equation (j2.30p as 

D{X;n,t) = det[I + B{i -V)A]. (2.34) 
It follows from Eqs. (l232l) and (f2:33]l that 

Bii-V)A=^l °), (2.35) 

where 

k = B21A12 + B22A22 - ^22V2lil2. (2.36) 

The expression on the right hand side of (j2.36p is evaluated with the help of identity 

—'n{k-q)ri{q-p)=r]{k-p) (2.37) 
-■K 27r 

which gives for the integral kernel of the operator K in Eq. (j2.36p 

K{k,p) = [j* z{k)z^^ (p) + jz^^ {k)z{p)]r]{k — p) 

. r dq 2, 



+ ^^*z-'{k)z~'{p) J -L^{k-qUq-p)z'{q). (2.38) 
The identity 

r]{k - q)ri{q - p) = r]{k - p) [ri{k - q) - ri{p - q)] (2.39) 
allows us to represent the last term in ()2.38p as 

- qHq-p)z\q) = {^i^^k) + z\p)] + f{k) - /(p)| v{k-p), (2.40) 
where 

Substituting Eq. (j2.40p into (j2.38p we arrive at the expressions (jl.7p - (jl.l2p for the 
integral kernel of the operator (|2.36p . Combining this result with Eqs. (j2.34p and 
(1235]) we arrive at Eq. (fTBD . 



3 Conclusion 

Using Wick's theorem we derived a Fredholm determinant representation, Eqs. ()1.6p - 
(jl.l2p . of the time-dependent correlation function (jl.ip of the Jordan- Wigner op- 
erator ()1.2p . Generalization to other operators of the Jordan- Wigner type, such 
as CjOj, and QjOj, is straightforward within this approach. The obtained rep- 
resentation can be used to investigate the charge dynamics in the antiferromag- 
netic [21 El m [5], spin-incoherent [HIEliaiaESlETlEElESlEOlEIlES], and spin- 
polarized [H El [71 [33] ID gases. It would also be interesting to generalize the 
calculations given here to the finite-temperature case, in particular, in view of the 
current interest to the anyon problem [141 [IS] . Another possible use of the deter- 
minant representation is to investigate the applicability range of the low-energy 
effective field theory (bosonization) description of the functions (jl.ip . The latter 
description is known [34j to suffer from Stokes-like phenomena such as the loss of 
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the A-periodicity of the large n, t asymptotics of the function (jl.ip . A way to get a 
A-periodic asymptotics could be the use of the technique based on the asymptotic 
solution of the matrix Riemann-Hilbert problem. This technique, developed in 
Refs. [351 EH EZ] proved to be very helpful (see, for example, [281 EB] and references 
therein) for the analysis of the asymptotic behavior of the Fredholm determinants 
with a kernel of the type (jl.Sp . 
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